Inertial Range Scaling in Rotations of Long Rods in Turbulence 
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We measure the rotational statistics of neutrally buoyant rods with lengths 2.8 < l/rj < 72.9 in turbulence. 
For particles with length in the inertial range, we derive a scaling relationship for the mean square rotation rate, 
{piPi) oc l~*'^ and show that measurements approach this scaling. Deviations from the proposed scaling are 
explained as the effect of dissipation range scales. The correlation time of the Lagrangian autocorrelation of 
rod rotation rate scales as the turn over time of eddies of the size of the rod. Measuring rotational dynamics of 
single long rods provides a new way to access the spatial structure of the flow at different length scales. 



The dynamics of particulate material in fluid flows is im- 
portant in a broad range of problems in nature [lT| and in- 
dustry 12(1 . Particles that are small and neutrally buoyant ex- 
hibit the same translational dynamics as fluid elements and are 
considered Lagrangian tracers. However, in many cases the 
particulate matter has different dynamics from fluid elements 
because it is not density matched with the fluid or has sizes 
larger than the smallest length scales of the flow. The mo- 
tion of density mismatched spherical particles in fluid flows 
has received a lot of attention 13|-|g] . Recently, the motion of 
neutrally buoyant, large, spherical particles in turbulent flows 
have become accessible in experiments iIt MioII and numerical 
simulations Illli[l2ll . and this has provided a novel method by 
which the spatial structure of turbulent flows can be probed 
with single particle measurements. The dynamics of spheres 
that are both large and density mismatched is more complex 
and has been addressed by some very recent experiments uM . 

Anisotropic particles have very different translational and 
rotational dynamics from spherical particles. Measurements 
and simulations of small rod-like particles show preferential 
alignment with the velocity gradients of the flow 11141 [1511 and 
this alignment suppresses the measured rotation rate in turbu- 
lent flows llm UTIl . A wide range of experimental and numer- 
ical studies have explored the dynamics of neutrally buoyant 
small rods and fibers in different flows [Il8i - l22ll . Only a few 
studies have focused on the dynamics of long rods in turbu- 
lence, where rod length s pan s over many times rj, the Kol- 
mogorov length scale jm 12311 . Using numerical simulations 
and slender body theory. Shin and Koch 111611 studied the trans- 
lational and rotational dynamics of long fibers in turbulent 
flow at Taylor Reynolds number up to Rx ~ 53.3. Among 
other things they show how the mean square rotation rate de- 
creases as rods become longer than the tracer limit and iden- 
tify the key role played by alignment of rods. 

It is known that the acceleration variance of large spheres 
scales with their diamete r app roximately as (a^) ^ dr"^!^ for 
d in the inertial range iItI-UOII. This result can be obtained by 
dimensional arguments simply assuming that the sphere ter- 
minates the cascade at its diameter Q. It can also be derived 
from the inertial range form of thepressure structure func- 
tion or the acceleration correlations llsl-lioll. This suggests that 
measuring the acceleration of large spheres provides access to 
the statistics of the turbulence at scales equal to the size of 



the spheres using only single particle measurements. Recent 
experimental work by Volk et al llioll finds that the dr'^l'^ scal- 
ing is not exact and proposes an improved scaling of d^"^^ 
by including intermittency in the pressure structure function. 
The fact that the accelerations of large spheres provides infor- 
mation about inertial range scaling suggests that similar infor- 
mation should be available in the rotations of long rods. 

In contrast to the point-particle models used for small in- 
ertial particles, the equations of motion for large spheres are 
largely unknown. Faxen corrections can be used to extend 
point particle models to describe large spheres |j24|, but these 
models have difficulties when particles are much larger than 

For rods, the analytical work can be done more rigorously 
to connect the motion of rods to the fluid motion. Olson and 
Kerekes ||23I1 have introduced a model to describe the rota- 
tional velocity of fibers. This model is based on the assump- 
tion that the fibers are infinitely thin and composed of a num- 
ber of sections smaller than 77 which are hydrodynamically 
independent. For a neutrally buoyant fiber of length I, the ro- 
tation rate is 
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where p is the orientation unit vector of the fiber and u is the 
turbulent velocity at points along the fiber If the orientation 
of a rod is uncorrelated with the velocity field in Eq. [T] then 
the mean square rotation rate of randomly oriented long rods 
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where u'^ is the rms velocity of the fluid flow and Rnn{s) is 
the fluid transverse velocity correlation function at separation 
distance of r. Shin and Koch lllql show that Eq.[2]is in good 
agreement with their simulations for the case of randomly ori- 
ented rods. 



In this paper, we introduce a scaling for the mean square 
rotation rate of rods with lengths in the Kolmogorov inertial 
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TABLE I: Table of flow parameters: Rx, Taylor Reynolds number; u, rms velocity of the flow; L, energy input length scale; (e), energy 
dissipation rate; L, energy input length scale; 77, Kolmogorov length scale; r^, Kolmogorov time scale, v is the fluid kinematic viscosity and 
isl.75xlO~'^m^/s. 



range. This scaling can be obtained either from dimensional 
arguments or from Eq. |2] The rotation rate has dimensions of 
inverse time, so for tracer rods it scales with the Kolmogorov 
time scale, (piPi) ^ t,,^^. Assuming that long rods are only 
rotated and aligned by eddies close to their size, the mean 
square rotation rate for rods at length scale / will scale like 
Ti~'^, where r/ is the time scale of eddies of size /. In the 
inertial range the time scale t; can be defined as r/ = l/ui = 
l/{l{e)Y^^), where ui is the velocity at length I, and (e) is 
the mean energy dissipation rate. This dimensional argument 
gives (jpiPi) '~ l^^/^ for / in the inertial range. 

The same l~'^/'^ scaling can be obtained using Eq.|2] with 
the additional benefit that the coefficient of the scaling law 
can be determined in terms of the coefficient of the second or- 
der structure function. In homogeneous turbulence the trans- 
verse correlation function is given in terms of the transverse 
velocity structure function, u^RMNir) = ii? — l/2D]^N{r). 
In the inertial range of isotropic turbulence, DNN{r) = 
(4/3)C2( (e)r)^^ /'^\ where C2 is an approximately universal 
constant ll25ll26ll . Using this form in Eq.|2] the mean squared 
rotation rate for randomly oriented rods in the inertial range is 
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where ly is the kinematic viscosity. Here we use the fact that 
only scales near the length of the rod contribute to the rotation 
rate in Eq. |2] so that in the high Reynolds number limit the 
inertial range form of the structure function can be used for 
all r. 

We have performed a series of experiments to measure the 
rotation of rods in 3D turbulence for rod lengths that ex- 
tend from the dissipation range well into the inertial range 
to explore whether an /^^/■^ scaling range exists. We mea- 
sure neutrally buoyant rod like particles with lengths ranging 
from 2.877 up to 72.977 in a turbulent flow between oscillating 
grids 1I27I1 . The rods are nylon thread with diameter of 0.2 
mm and are cut to different lengths (/ =1, 3, 6.8, 15.2 mm). 
All particles are dyed fluorescent for better detection. The 
rotational dynamics of rods are measured using stereoscopic 
images from four high speed cameras IUTIi . Each camera that 
detects a rod defines a plane in which the rod should exist. 
We require at least three cameras to detect each single rod 
and the orientation of the rod is the intersection of the planes 
these cameras define. The detection volume is illumined with 
4 laser beams. This has nearly removed a limitation in earlier 
experiments IUTII where the probability of detecting a particle 



depended on the orientation of the particle with respect to the 
laser beam. 

The rotation rate vector of rods, p, is measured from 
quadratic fits to the measured orientations along trajectories 
versus time. The measurements are at two different Taylor 
Reynolds numbers {R\ =150, 210). We have done parallel 
experiments with tracer particles to measure the turbulence 
parameters. From measured tracer velocities, we extract the 
second and third order longitudinal structure functions and 
obtain the energy dissipation rate, (e), from Kolmogorov's 4/5 
law. We obtain the flow parameters shown in Table Ufrom the 
measured energy dissipation rate and rms flow velocity. 

In these experiments the maximum detection volume is 160 
cm'^, while the effective detection volume is different for each 
rod length and is the smallest (85 cm'^) for the longest rods at 
15.2 mm. The effective detection volume is smaller to ensure 
that the entire rod is in the illuminated detection volume so the 
position and orientation of rods are measured more accurately 
based on the full length of rods. The number density of rods is 
small so particle-paiticle interaction is negligible. The particle 
concentration is 0.025 cm~'^ for 1 mm rods, and 0.0075 cm~'^ 
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FIG. 1 : (Color online) The PDF of rotation rate squared for different 
rod lengths. The lengths of rods are l/ri= 2.8 (green crosses), 4.9 
(black filled circles), 8.5 (right triangles), 14.5 (open squares), 19.1 
(red open circles), 32.8 (blue diamond), 42.2 (brown left triangles) 
and 72.9 (purple asterisk) and tracers from simulation (solid gray 
line). The results are reported from two experiments at Rx — 150 
and 210. The simulations ! 1711 are for tracer rods at i?A=180. 



for the longest rods at 15.2 mm. The uncertainty in measuring 
the center of rods is 60/xm for the 1 mm rods and 180 /^jm for 
15.2 mm rods. This uncertainty is determined based on the 
stereomatching accuracy. The accuracy of measuring the ori- 
entation for 1 mm particles is 0.01 rad and increases with rod 
length. The uncertainty in the orientation of the rods is deter- 
mined from the residual of the intersection of planes defined 
by multiple cameras. 

Figure [T] shows the probability distribution function (PDF) 
of the rotation rate squared, piPi, normalized by the mean for 
different rod lengths. The PDF shows only a weak depen- 
dence on the rod length and all rod lengths show rare events 
with large rotation rates. The probability of rare events is 
somewhat smaller for long rods {l/r] > 20) than for tracer 
rods {l/rj < 7); however, this difference is slightly larger than 
the measurement uncertainty due to the smaller number of 
samples for long rods. The error bars represent the random 
statistical error and the systematic error in measuring the ro- 
tation rates. Qualitatively, the rotation rate PDF depends on 
rod length in the same way that the acceleration PDF for large 
spheres depends on diameter f\M, with only a small narrow- 
ing of the tails for large particles. 

The mean square rotation rate for different rod lengths is 
shown in Fig.|2a). Our experimental measurements show that 
by increasing the length of rods, the mean square rotation rate 
decreases. This is expected since longer rods should begin to 
filter out the contributions from some eddies smaller than their 
length. These results agree with an earlier simulation [16J of 
long fibers in turbulent flow at R\ =53.3. The experimental 
measurements are slightly larger, but the difference is only a 
little larger than the measurement uncertainty. 

For a more detailed understanding of the dependence of the 
rotation rate on the length of rods, we compare these results 
with the rotation rate of randomly oriented rods. Assuming 
the orientation of rods is uncorrected with the fluid veloc- 
ity, Eq. 12] can be used to calculate the mean square rotation 
rate from measured second order structure functions. Fig- 
ure 12 a) shows the mean square rotation rate for randomly 
oriented rods predicted by this model using our experimen- 
tal measurements of D^^n at i?A=150, 210. The mean square 
rotation rate of short rods is much smaller than randomly ori- 
ented rods of the same length {l/rj). However, this difference 
decreases as the length of the rod is increased. Previous stud- 
ies of tracer rods 11141 Im UTI l28l 12911 have shown that as rods 
are carried by the flow their orientation becomes correlated 
with the directions defined by the velocity gradient tensor of 
the flow. This alignment results in suppression of the rota- 
tion rate of short rods compared to randomly oriented rods 
lUTIl . The smaller differences between the measured rotation 
rates of long rods and randomly oriented rods suggests that the 
alignment is weaker between long rods and the velocity gra- 
dients responsible for rotating the long rods. In the previous 
simulations of long fibers lUal the correlation of long fibers 
with the flow is much weaker than what we see in the ex- 
periment, possibly as a result of their much smaller Reynolds 
number. 

Figure lltb) shows the inertial range scaling law from Eq.|3] 
with 6*2 — 2.0 and compares it with the experimentally mea- 
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FIG. 2: (Color online) Mean square rotation rate as a function of rod 
length, a) Comparison of the experimental data at Rx =150 (red open 
circles) and 210 (red squares) with the model of randomly oriented 
rods at i?A=150 (puiple diamond) and R\= 210 (black triangles), b) 
Comparison of the mean square rotation rate of rods (red open circles 
and squares) and the model for randomly oriented rods (diamond and 
triangles) with an Z^''' "^ inertial range scaling law (solid green line). 



sured rotation rates and the prediction of Eq. |2]for randomly 
oriented rods using the measured velocity structure function. 
Both the measured rotation rates and the prediction of Eq. |2] 
approach an /^^/'^ scaling for large /. The experimental data 
has a smaller coefficient compared to randomly oriented rods 
as expected due to alignment effects. Within error bars of the 
experimental data for I > SO??, one could also fit the data with 
a different exponent slightly steeper than —4/3. In this same 
range , the prediction of Eq. |2]is also steeper than l^^^'^. The 
cause of the steeper scaling can be found in the fact that the 
prediction of Eq. |2] overshoots the power law scaling in the 
range 20?7 < I < 5Qrj. This overshoot occurs in the range of 
scales slightly larger than the dissipative range because scales 
smaller than the length of the rod contribute to the rotation 
rate. For / > 5O77, well into the inertial range, the prediction 
converges with the scaling law from Eq. |3] because here the 
contributions from the dissipation range are becoming negli- 
gible. The overshoot also appears to exist in the experimental 



rotation rate data, although the effect is of the same size as the 
measurement uncertainties. 

The sUghtly steeper scaling of the rod rotation rate is sim- 
ilar to the effect observed in Ref. |ldl for the accelerations 
of spheres with diameters in the inertial range. They argue 
that intermittency effects are responsible for the difference 
between their measured scaling and the prediction of dimen- 
sional analysis. It is possible that intermittency effects also 
play a role for rods. However, the availability of a solid the- 
oretical foundation for calculating rotation rates of long rods 
suggests another possible explanation due to effects of small 
scales. We have used the Batchelor parameterization IBOl 13111 
of the structure function at very high Reynolds number in 
Eq. ID and find an overshoot which leads to steeper scaling 
for 3O77 < I < IOO77 and then an agreement with /^^/■^ for 
I > lOOrj. This parameterization has no intermittency effects. 
It is possible that spheres also have an overshoot in which par- 
ticles only a little larger than the dissipation scale have accel- 
erations larger than the inertial range power law due to the 
contributions from dissipation range scales. 

To measure the mean square rotation_rate, we use the ex- 
trapolation method developed in Ref. ||7|] to remove the de- 
pendence on the smoothing used to measure derivatives. The 
extrapolation to zero smoothing length is known to overesti- 
mate the actual rotation rate, and we performed a simulation 
of our measurement system to estimate the size of this over- 
estimate. The translational motion of tracer rods matches that 
of fluid particles, so we are able to use a database of previ- 
ously simulated Lagrangian trajectories [|32|| to integrate Jef- 
fery's equation [33] and obtain the orientation of tracer rods 
along Lagrangian trajectories. We have used these trajectories 
to create simulated experimental images of rods and analyzed 
these images with the same analysis codes used for the ex- 
perimental data. The extrapolated mean square rotation rate 
of simulated rods is found to be larger than direct measure- 
ments of the DNS rods by 6%. We have corrected for this by 
shifting our mean square rotation rate measurements down by 
6%. The error bars on the experimental measurement of the 
mean square rotation rate in Fig. |2]represent both statistical 
uncertainty and the systematic error due to the overestimation 
correction. 

The small difference between two Reynolds numbers in the 
model for the mean square rotation rate of randomly oriented 
rods in Fig. |2]is due to uncertainties in measuring the sec- 
ond order transverse velocity structure function. As there are 
very few tracer particle pairs measured with very small sep- 
arations (r < 15rj), we used analytical forms of the struc- 
ture function there. For r < ?/, the dissipation range is 
DNN{r) = j^{{e) / i')r'^ . For r] < r < 15i], we use a fit 
of the Batchelor parametrization OOl [3l1l . This method pro- 
vides a smooth transition from small scales to the region of 
experimentally resolved velocity structure function. 

Figure |3la) shows the Lagrangian autocorrelation of rota- 
tion rate measured for different rod lengths. Our measure- 
ments show that the correlation time of the rotation rate de- 
pends on the length of rods and increases with rod length. 
Similar autocorrelation functions were obtained in simula- 
tions 1,16,1 at _Ra=39.9. We expect that if the rods are rotat- 
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FIG. 3: (Color online) Lagrangian autocorrelation of rotation rate 
for different rod length. The lengths of rods are l/ri= 2.8 (green 
crosses), 4.9 (black filled circles), 8.5 (right triangles), 14.5 ( open 
squares), 19.1 (red open circles), 32.8 (blue diamond), 42.2 (brown 
left triangles) and 72.9 (purple asterisk), a) Time is normalized by 
the Kolmogorov time, b) Time is normalized by the turn over time 
of eddies with size equal to the length of the rods, r;. The symbols 
are displayed at every other data point. 



ing due to eddies of their size, /, then the decay time for the 
correlation of rotation rate should scale as the turn over time 
at the length of the rods. In Fig. |3lb) the horizontal axis (f) 
is normalized by t;, time scale of eddies with length scale 
/. After this normalization, the Lagrangian autocorrelation 
of rotation rate for all rod lengths collapse on a single curve 
within measurement uncertainty. The time scales for short 
tracer rods (l/ij < 5) is the Kolmogorov time-scale (r,;) and 
the time scales for longer rods are measured from longitudinal 
second order velocity structure function of the fluid particles 
in = -^l/Sui = -^1/^Dll{1))- Measuring the autocor- 
relation function for long rods (/ = 42. 4?] and 72.977) is dif- 
ficult as the effective detection volume is small compared to 
the length of these rods so the trajectories are not long enough 
to measure long time autocorrelations. The collapse seen in 
Fig. [3jb) provides additional evidence, beyond that seen in 
the mean square rotation rate data, that rod rotations are con- 
trolled by eddies with size near the rod length. 



Long rods provide a promising path for studying the dy- 
namics of large particles in turbulence. For rods, slender body 
theory can be used to connect particle motion with the fluid 
flow even for particles much larger than the Kolmogorov scale 
while the analytical results for spheres are only available for 
small deviations from the small particle limit. We have pre- 
sented measurements of the rotations of rods with lengths ex- 
tending well into the inertial range, and we find that the mean 
square rotation rate in the inertial range approaches the l~^/^ 
scaling that is expected from inertial range scaling of the ve- 
locity structure functions. The PDF of rotation rate shows 
only a weak dependence on rod length. We find that rods 
develop preferential alignment so that their rotation rates are 
significantly smaller than that predicted for randomly oriented 
rods. This alignment depends on rod length as rods in the iner- 
tial range show a smaller effect of alignment than tracer rods. 
The Lagrangian autocorrelation time of the rotation rate de- 
pends on the length of rods and scales with the eddy turn over 



time at a scale equal to the rod length. 

Shin and Koch flq] provided a groundbreaking simulation 
data set on this problem, but were limited to R\ < 53.3 where 
there is essentially no inertial range. Future high Reynolds 
number DNS using their method of simulating long fibers of- 
fers the possibility to study the motion of long rods while also 
having access to the full velocity field around the rods. Ex- 
perimental tracking of long rods in turbulence allows access 
to the dynamics of turbulent scales at the length of the par- 
ticle from single particle measurements, and has potential to 
provide valuable information about Lagrangian dynamics as a 
function of scale in complex turbulent flows. 

We acknowledge support from NSF grant DMR- 1208990, 
and COST Actions MP0806 and FP1005. We thank Stefan 
Kramel for assistance in the experimental work and Enrico 
Calzavarani, Federico Toschi, Nicholas T. Ouellette and Rui 
Ni for stimulating discussions. 



[1 

[2 

[3 

[4 
[5 

[6 

[7 

[8 

[9 

[10 

[11 
[12; 

[13 



[14 
[15 

[16 
[17 



Bowen, J. D., Stolzenbach, K. D., and Chisholm, S. W. Limnol- 
ogy and Oceanography 38(1), 36-51 (1993). 
Lundell, R, Soderberg, D., and Alfredsson, H. Annii. Rev. Fluid 
Mech. 43, 195-217(2011). 

Falkovich, G., Fouxon, A., and Stepanov, M. G. Nature 
419(6903), 151-154 (2002). 

Shaw, R. A. Annu. Rev. Fluid Mech. 35, 183-227 (2003). 
Toschi, F. and Bodenschatz, E. Annu. Rev. Fluid Mech. 41, 
375-404 (2009). 

Meiburg, E. and Kneller, B. Annu. Rev. Fluid Mech. 42, 135- 
156 (2010). 

Voth, G. A., La Porta, A., Crawford, A. M., Alexander, J., and 
Bodenschatz, E. J. Fluid Mech. 469, 121-160 (2002). 
Qureshi, N. M., Bourgoin, M., Baudet, C, Cartellier, A., and 
Gagne, Y. Phys. Rev Lett. 99, 184502 (2007). 
Brown, R. D., Warhaft, Z., and Voth, G. A. Phys. Rev Lett. 103, 
194501 (2009). 

Volk, R., Calzavarini, E., Leveque, E., and Pinton, J.-F. / Fluid 
Mech. 66», 222,-21,5 {2011). 

Homann, H. and Bee, J. J. Fluid Mech 651, 81-91 (2010). 
Calzavarini, E., Volk, R., Leveque, E., Pinton, J. F, and Toschi, 
F. Physica D-Nonlinear Phenomena 241(3, SI), 237-244 
(2012).' 

Prakash, V. N., Tagawa, Y., Calzavarini, E., Mercado, J. M., 
Toschi, F, Lohse, D., and Sun, C. New J. Phys. 14 105017 
(2012). 

Pumir, A. and Wilkinson, M. New J. Phys. 13 093030 (201 1). 
Parsa, S., Guasto, J. S., Kishore, M., Ouellette, N. T, GoUub, 
J. P, and Voth, G. A. Phys. Fluids 23(4) 043302 (201 1). 
Shin, M. and Koch, D. J. Fluid Mech. 540, 143-173 (2005). 
Parsa, S., Calzavarini, E., Toschi, F, and Voth, G. A. Phys Rev. 



Lett 109(13) (2012). 
[18] Parsheh, M., Brown, M., and Aidun, C. J. Fluid Mech. 545, 

245-269 (2005). 
[19] Holm, R. and Soderberg, D. Rheologica Acta 46(5), 721-729 

(2007). 
[20] Lin, J., Zhang, S., and Olson, J. A. Engineering Computations 

24(1-2), 52-76(2007). 
[21] Mortensen, P. H., Andersson, H. I., Gillissen, J. J. J., and 

Boersma, B. J. Phys. Fluids 20(9) 093302 (2008). 
[22] Marchioli, C, Fantoni, M., and Soldati, A. Physics of Fluids 

22(3) 033301 (2010). 
[23] Olson, J. and Kerekes, R. J. Fluid Mech. 377, 47-64 (1998). 
[24] Calzavarini, E., Volk, R., Bourgoin, M., Leveque, E., Pinton, 

J. F, and Toschi, F J. Fluid Mech 630, 179-189 (2009). 
[25] Pope, S. B. Turbulent Flows. Cambridge University Press, 

(Cambridge), (2000). 
[26] Sreenivasan, K. R. Phys. Fluids 7(11), 2778-2784 (1995). 
[27] Blum, D. B., Kunwar, S. B., Johnson, J., and Voth, G. A. Phys. 

F/Mrf.? 22(1) 015107 (2010). 
[28] Girimaji, S. S. and Pope, S. B. J. Fluid Mech. 220, 427^58 

(1990). 
[29] Luthi, B., Tsinober, A., and Kinzelbach, W. J. Fluid Mech. 528, 

87-118(2005). 
[30] Batchelor, G. K. Proc. of the Cambridge Philosophical Society 

47(2), 359-374(1951). ' 
[31] Grossmann, S. Phys Rev E 51(6, Partb), 6275-6277 (1995). 
[32] Benzi, R., Biferale, L., Calzavarini, E., Lohse, D., and Toschi, 

F Phys Rev E 80 066318 (2009). 
[33] Jeffery, G. B. Proc. R. Soc. Lond. A 102, 161-179 (1922). 



